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Abstract 

For a class of polyhedrons denoted K n (r, e), we construct a bijective contin¬ 
uous area preserving map from K n (r, e) to the sphere § 2 (r), together with 
its inverse. Then we investigate for which polyhedrons K n (r',£) the area 
preserving map can be used for constructing a bijective continuous volume 
preserving map from K n (r',e) to the ball § 2 (r). These maps can be further 
used in constructing uniform and rehnable grids on the sphere and on the 
ball, starting from uniform and rehnable grids of the polyhedrons K n (r, e) 
and Kjgr'je), respectively. In particular, we show that HEALPix grids can 
be obtained by mappings polyhedrons K n (r, e) onto the sphere. 


1. Introduction 

A uniform grid on a two(three)-dimensional domain D is a grid all of 
whose cells have the same area(volume). This is required in statistical appli¬ 
cations and in construction of wavelet bases of the space L 2 (D), when one 
wishes to use the standard inner product and 2-norm, instead of a weighted 
norm dependent on the grid [§j]. A refinement process is needed for multires- 
olution analysis or for multigrid methods, when a grid is not fine enough to 
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solve a problem accurately. A refinement of a grid is called uniform when 
each cell is divided into a given number of smaller cells having the same 
measure. To be efficient in practice, a refinement procedure should also be 
a simple one. In many applications, especially in geosciences, one requires 
simple, uniform and reffiiable (hierarchical) grids on the 2D-sphere or on the 
3D-ball (i.e. the solid sphere). One simple method to construct grids on the 
2D-sphere is to transfer existing planar grids, while for the 3D-ball a simple 
way would be to transfer polyhedral grids. 

Partitions of the 2D-sphere § 2 into regions of small diameter and equal 
area has already been constructed by Alexander [l], In [g| Leopardi derives a 
recursive zonal equal area partitioning algorithm for the unit sphere E> d em¬ 
bedded in M d+1 . The partition for the particular case d = 2 consists of polar 
cups and rectilinear regions that are arranged in zonal collars. Besides the 
fact that that the regions have different shapes, his partition is not suitable 
for applications where one must avoid that vertices of spherical rectangles lie 
on edges of neighbor rectangles. 

In astronomy, the most used construction of equal area partitions of S 2 , 
is the HEALPix grid [ 2 ], providing a hierarchical equal area iso-latitude 
pixelization. Other constructions are the truncated icosahedron-method of 


141, the 


10 


Snyder [13], the small circle subdivision method introduced in 
icosahedron-based method by Tegmark (Toj , see also 0- Section 1 in 
contains a larger list of uniform spherical grids, together with their proper¬ 
ties. A complete description of all known spherical projections from a sphere 
or parts of a sphere to the plane, used in cartography, is realized in [ 3 , 15]. We 
should mention that most of the existing constructions of spherical hierar¬ 
chical (i.e. refinable ) grids do not provide an equal area partition. However, 
hi ]91, 10, [HI we have already constructed some area preserving maps onto 
the sphere, using some Lambert azimuthal equal area projections. By trans¬ 
porting onto the sphere uniform and reffiiable planar grids, we could obtain 
uniform and refinable spherical grids. Also, in ( 4 ] we have constructed an 
octahedral equal area partition of the sphere without make use of Lambert 
projection, and we have studied some properties of the corresponding config¬ 
urations of points. Regarding the equal volume partition of the 3D-ball, to 
our knowledge there was no uniform and reffiiable grid. This was the motiva¬ 
tion for the work in 12], where we have constructed a volume preserving map 
from the 3D-cube to the 3D-ball. This allowed us not only the construction 
of uniform and refinable grids on the 3D-ball, but also a uniform sampling 
in the space of 3D rotations with, with applications in texture analysis. 
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In this paper we consider a class of convex polyhedrons K n (r, e) and first 
we construct an area preserving map from K n (r, e) to the 2D-sphere of radius 
r. Then, using this map, we construct a volume preserving projection from 
the interior of the polyhedrons to the 3D-ball of radius r. The spherical grids 
in [ 2 ] and Q can be obtained by mapping grids on K n (r, e) to the 2D-sphere. 

The paper is structured as follows. In Section [ 2 ] we describe the class of 
polyhedrons K n (r, e) which will be projected onto the sphere and we give 
some formulas which will be useful in the next sections. In Section [3] we 
construct a continuous area preserving map from K ra (r, e) to the sphere, 
together with its inverse. In Section 13.51 we show how our maps can be 
useful for the construction of uniform and rehnable grids on the sphere. In 
particular, we show that the HEALPix grids [ 2 J can be obtained as continuous 
images of polyhedrons K n (r, e) . Finally, in Section [4] we construct a volume 
preserving map from the interior of the polyhedrons K n (r, e) to the 3D-ball 
of radius r, using the new constructed area preserving map. 

2. Preliminaries 

Let r > 0 and consider the sphere § 2 (r) = § 2 of radius r centered at the 
origin O, of parametric equations 

x = r cos 0 sin <^, 

y = r sin 6 siny?, (1) 

z = r cos 99 , 

where p E [0,7r] is the colatitude and 6 G [0, 27t) is the longitude. A simple 
calculation shows that the area element of the sphere is 

dS = r 2 sin p d6 dp. (2) 

Let e G [0,1) be fixed. We intersect the sphere S 2 with the planes z = 
±£T, and we denote with = § + the spherical cap situated above the plane 
z = er, with §“ = S - the spherical cap situated below the plane z = —er, 
and with E er = E the equatorial belt between the planes z = ±er. Their 
areas are 

*4.(§ + ) = *4.(§“) = 27t( 1 — e)r 2 , *4.(E) = 47rer 2 , 

and since § + U §” U E = § 2 and these portions are pairwise disjoint, one has 

^(§+) + A(S~) + A(E) = ^(§ 2 ) = 4vrr 2 . 
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Figure 1: The polyhedron K„(r, e) which will be projected on the sphere. Here n= 6. 

For a fixed integer n, n > 3, let K n = K n (r, e) be a polyhedron formed 
by a regular prism B n = B n (r, e) of height 2 er and two congruent pyramids 
P+(r, e) = P+ and P“(r, e) = P“ of lateral area 27 t(1 — e)r 2 , such that the 
two bases of the prism B n coincide with the bases of the two pyramids (see 
Figured]). Let V n denotes the regular polygon which is the base of P+ and the 
upper base of B n . We place K ra such that B n is centered at O and symmetric 
with respect to the plane OXY and such that one of its vertical edge is 
situated in the plane OXZ. Further, we divide the space R 3 onto 2 n zones 
if, i = 0, ..., n — 1, as follows: 

If = {( x,y,z ) ei 3 y>0, 0 < y < xtan^} , 

= {Of 2/,-) e M 3 ,^ < 0, 0 <y < Xtan^} , 

and further, the other zones if are obtained by rotating if with the angle 
a.i = — as 

6 n 

if = {Ki ■ (x, y> z ) T , (x, y , z) e if }, 

where IZ t is the 3 1) rotation matrix around OZ, 

( cos cc, — sin a* 0 
sin a, cos a.i 0 

0 0 1 

Thus, each face of the pyramids Pj will be situated in one of the domains 
if (see Figure [2]). 
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Figure 2: The portions Iq and I 0 . 


In Section [3] we will construct a map T n , r = T n '■ § 2 —* K n which preserves 
areas, in the sense that 

A{D) = A{T n {D)), for all D C § 2 , (3) 

where A(D) denotes the area of a domain D. For an arbitrary point (x, y, z) 6 
S 2 we denote 

{X,Y,Z) = T n (x,y,z)eK n . (4) 

The restrictions of T n to § + , S _ , E will be denoted by 7^ + , T~ and 
T^, respectively. More precisely, we will deduce the formulas for the area 
preserving maps 7^ + : § + —> P+, TT :§”—)• and : E + —> B n . 

We introduce the following notations: 

R n = the radius of the circle circumscribed to V n , 
r n = the radius of the circle inscribed in V n , 
b n = the altitude of the pyramid P^, 

£ n = the edge of the polygon V n , 
a n = the slant height of the pyramid P^. 


We impose that the lateral area of B n equals the area of the equatorial belt 
E, therefore 



( 5 ) 


n 

On the other hand, since £ n = 2 R n sin -, we deduce that 

_ vrr 
n n sin - 


( 6 ) 


n 
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We also have 


r n = R n cos . 

n 


The area of a face of the pyramid P+ is 



(7) 

( 8 ) 


and by imposing that *4.(P+) = *4.(§ + ), we find that a n = 2(1 — e)r and 


An — 2^(1 ~ £ )r 2 — £ n (l — £ )r . 


n 


( 9 ) 


Using the equality a„ = rf + b„, we obtain 


h n — r[ 4(1 — j cot 2 —^ 
' n z n) 


( 10 ) 


3. Construction of the area preserving map 'T n : § 2 —> ¥L n and its 
inverse 

3.1. Construction of the map Tff 

For h G (1 — e, 1), we denote by E> r h the spherical cap situated above the 
plane of equation z = rh. A simple calculation shows that 

-4(§r-h) = 27 t(1 — h)r 2 . (11) 


Now, we calculate H > 0 such that the portion P f(H) of the pyramid P+ 
situated above the plane z = H has the same area A(E> r h). For the small 
pyramid P f(H), let R n be the radius of the circumscribed circle of its base. 
Then, one has 

R n _ K - H + er 

D T ’ 

ll n u n 

and the area A(H) of a face of P f(H) is therefore 


A n (H) 


A 



H 


er 


2 


27t(1 — e)r 2 f H — er\ 2 

n V K ) 
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Imposing now A(E> r h) = nA(H), we obtain that 


2 


27r(l — h)r 2 = 27r(l — e)r 2 



H — er 
bn 


whence after some calculations we obtain 


H 


er + b 


n 



l -h 
1 — e 


In conclusion, for z > er we can define Z from (J4J) as 


Z = er + b n 



1 — z/r 


1 — e 


( 12 ) 


meaning that the parallel circles of the sphere map onto polygons obtained 
as intersections of P+ with planes parallel to XOY. 

For the formulas for X and Y we proceed as follows. We focus on the 
portion Xq of § 2 , situated in Iq, and the face of the pyramid P+ situated 
in Iq will be denoted with F 0 + . We denote A{\/ 1 — e 2 r , 0, er ), C(0,0,r) and 
we consider the vertical plane of equation y = xtancx, with a G (0, 27r/n) 

(see Figure [3j left). We denote by MC its intersection with Xq . More pre¬ 
cisely, M(y/1 — e 2 r cos a, y/1 — e 2 r sin ck, er). The area of the spherical do¬ 
main AMC equals a(l — e)r 2 . Now we intersect the face Fq of the pyramid 
with the vertical plane of equation y = xtan/1 and we denote by M'(X', Y') 
its intersection with the edge A'B', where A'(R n , 0, er) and B'(R n cos 2tt/ n, R n sin 27r/n, er) 
(see Figure [31 right). 

For finding the coordinates X 1 , Y' we use the fact that M' is the inter¬ 
section of the line A'B' with the plane of equation Y' = X' tan/3, so solving 


the system 



we obtain 


1 + tan /3 tan ^ ’ 1 + tan /3 tan ^ 


Rn y> - Rn taI1 ^ 


(13) 


Further, 
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Figure 3: The portion AMC of the spherical cap § + and its image triangle A'M'C' = 
T n (AMC) on the pyramid P+. 


and therefore, for the area of the planar triangle A'C'M' we find 


A(A'C'M') 


A'M' ■ a n 
2 


27r(l — £)r 2 ta,n/3 
n sin ^ (l + tan /3 tan 


If we impose that this area equals the area of the spherical domain AMC , 
which is a(l — e)r 2 , we finally find that 


tan /3 


• r 

an sin — 

_ n _ 

2 ( 7 r — an sin 2 


(14) 


This is in fact the relation between a and (3, such that the slice ACM of J r 0 1 ' 
has the same area with the triangular slice A'C'M' of the face P 0 + . 

Consider now V = (m, q,p) G 3Fq and denote N = ( u,s,t ) := T n (V). The 
point V is the intersection of 3Fq with the planes 


(Pi) : z = p and (P 2 ) : y = —x, 

m 


where — = tana. 
m 


The point N will be the intersection of T n (Pi), T n {P 2 ) and the face F 0 + , 
therefore 


T n (Pl) 

r n (p 2 ) 

^ 0 + 


t — er 


1 — p/r 

1 — e 

n sin — • arctan — 


s = 


2{ji — n arctan 


— sin 2 

m n) 


U, 


7T 


u ■ b n + s ■ b n tan - + t ■ R n = (b n + £r)R n 
n 














Solving this system we find 


u 

s 

t 


R, 


R? 


1 — p/r 

1 — £ 


n . 2 7r q 

1-sin — arctan — 


7r 


n 


m 


1 — p/r n sin — arctan — 

r! n m 


1 — e 


2t r 


er + b n 1 - 


1 — p/r 
l — £ 


In conclusion, using (EJ) and with the notations in ([3]), for ( x,y,z ) G J r 0 f , we 
define the restriction of to /Fq as 


X 

Y 

Z 


r(r 


7r 


n sm - 


7r y 

sm — arctan — 


n 


r(r — z) 7r y 

-cos — arctan —, 

1 — £ n x 


£r 


1 - 


1 - z/r 

1 — £ 


X 


(15) 

(16) 
(17) 


with b n given in (fTOlh 

In general, for i — 0,..., n — 1, let J-\ be the portion of § 2 situated in 
if . The the restriction of Tf~ to /Ff will be calculated as 

(X,Y,Zf = 7ti-% + {nf ■{x,y,z) T ), for (x, y, z) G /Ff , he. (18) 


X 

Y 

Z 


r r 


1 — £ 


r r 


1 - 

£r + b„ 


7 T cos 7T -x sm on + y cos . 

---sulfa,; H —) ■ arctan-Hlfl) 

n sm - n x cos + y sm a, 

7 T sin a; . 7r —x sin a, + y cos a, , . 

- f + cos ai + - • arctan--—--- 

n sm - n x cos a, + y sm a. 


1 - 


1 — z/r 

1 — £ 


( 21 ) 
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In spherical coordinates (<p,9), the map 7p restricted to Eq writes as 


X 

Y 

Z 


7 r 


1 — e \ n sin 


a ■ M ■ <P 
— 9 sm — sm —, 


n 


2 a 71 ■ V 
- 9 cos — sm —, 

l — s n 2 


er + b n 1 


2 p 

■ sm — 


l-£ 2 


For the portion Ef the map is dehned by 


X 

Y 

Z 


1 — e \ n sin - 

' n 


7i sm op 


1 — £ \ n sin - 

' n 


7TCOSCp Tt(2i + 1) 

Sill- 

n 

7t(2 * + 1) 


+ cos 


n 


(6 - ) sin|, 

(9 - a^ ) sin^, 


£r + b n 1 - 


2 p 

-Sill — 

1 -£ 2 


If we evaluate the coefficients of the first fundamental order of this surface 


E' = pg 2 + (Yp 2 + (zp 2 = 




2(1-e) 

‘1 


4(1 — £) 2 + [e 


(2 i + 1)7r 


n 


cos 


F' = x' v x' e + rtf + z;z' = - 

G' = M 2 + (Yg ) 2 + (Z ') 2 = 


(2f + l)vr^ . 

9 -I sm</?, 


n 


2 r 2 . 2 V 

- sin —, 

l-£ 2 


then we find that E'G' — (E') 2 = r 4 sin 2 p, which is equal to EG — E 2 for the 
sphere (see (J2]) ). This property holds for all the portions Ef, in conclusion, 
a JEG — F 2 is invariant under 7p, whence 7p is an area preserving map. 

Remark 1. The map 7p : S _ — > P“ can be obtained by symmetry. More 
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precisely, the formulas are 


X 

Y 

Z 


r(r + z ) ( 7i 


1 — £ 


n sin ■ 


7 r y 

— sin — arctan — 
n x 


r(r + z) Ti y 

-cos — arctan —, 

1 — e n x 


er 


1 + z/r 


( 22 ) 

(23) 

(24) 


3.2. The inverse of the map Tff 

Let (X,Y,Z) e IP+. First, from (I2T]) we immediately deduce that 


z = r — r(l — e) 


Z -er\ 2 

b n ) 


Then, from (]T9]) and ([20]) we find that 


Y 

X 


7rsinQi 
nsin — 


+ cos (oti + ^) • arctan 


—x sin OLi+y cos a* 
x cos oti+y sinai 


7rcos aj 
n sin — 

n 


— sin (a* + ^) • arctan 


—x sin a.j -\-y cos a.j 5 
x cos oti -\-y sin oti 


(25) 


whence 


—x sm oti + y cos cxj 

arctan- 

x cos oti + y sin a* 


ti —X sin cxj + Y cos aj 

n sin - A" cos(aj + -) + Y sin(«j + -) 

n \ L n J v L n J 


Therefore, 


tan A cos aj + sm a* 

y = x ---, 

cos aj — tan A sm aj 

and then, taking into account that x 2 + y 2 + z 2 = r 2 , after simple calculations 
we find 


x = Vr 2 — z 2 ■ cos 
y = \/r 2 — z 2 ■ sin 
where z is given in (123]) . 


71 

—X sin aj + Y cos a* 

n sin - 

n 

X cos(aj + ^) + Y sin (a* + £) 

71 

—X sin aj + Y cos aj 

n sin - 

n 

A cos(aj + ^) + Y sin(aj + £) 



+ OLi 


(26) 

.(27) 
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Remark 2. For the inverse of T n , the formulas for x and y are as in (125]) 
and (1271) , respectively, while for z formula (1251) changes into 


z = r( 1 - e) 


— r. 


(28) 


Remark 3. Formulas (1261) and (1271) can also be obtained from the formulas 
for i = 0, similarly to (fl8l) . 

3.3. Construction of the maps (Tf) 1 and Tf 

The prism B n is regular with height 2 er and edge £ n given in flHJ). The 
face B 0 of B„ situated in the zone I 0 — Iq U If has the equation 

7T 7T 1 1 . , 

Y ■ sin —b X ■ cos — = R n cos —, (29) 

n n n 

or, equivalently, Y = ( R n — A") cot Then, in order to project the face B 0 
onto the plane OYZ , we make a translation with R n along the axis OX, 
followed by a rotation of angle — ^ around the axis OZ. Thus, the point 

{X, (Rn- X) cot -, Z) G B n 
n 

will be first translated to 


(X - Rn, (Rn - X) cot -, Z) G B n , 

n 

and then rotated around OZ it will be mapped onto the point 


P 0, 


R n - X 


sm - 

n 


,z) =p o, 


Y 


cos - 

n 


-,z 


As a translation followed by rotation, this map will be area preserving. Next 
we use the inverse Lambert cylindrical equal area projection (see i.e. |3j), 
the point P being further mapped onto 


\/r 2 — Z 2 cos-—, y/r 2 — Z 2 sin-, Z ) G E C S 2 . 


r cos ■ 


r cos 


(30) 


In conclusion, the map (Tf) 


e\-l . 


—>■ E is area preserving and maps the 


point (X, Y, Z) G B„ onto the point given in (j30]) . 
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In general, similarly to (I18p . the face of the prism B n situated in Jj = 
n u/f is mapped on the sphere by 


x 


= Vr 2 - Z 2 


cos 


y = \Jr 2 — Z 2 sin 
z = Z. 


Y cos aii — X sin a* 

r cos - 

n 

Y cos Oii — X sin a* 


r cos ■ 


+ ot, 


+ ot, 


For the direct application 7^ e : E —> B n , the calculations give 


(31) 

(32) 

(33) 


X 

Y 

Z 


„ . , 7T —X Sill OLi + y COS Oti 

R n cos a.i — r sin fa* H—) arctan-, 

n x cos oti + y sin a* 

. .7 r —x sin oti -\- y cos oti 

n ri sm oti + r cos (a* H—) arctan-, 

n x cos oti + y sin a* 

z. 


(34) 

(35) 

(36) 


3-4- The continuity ofT n 

Proposition 4. The map T n : § 2 —> K n is continuous. 


Proof. It is enough to restrict ourselves to the first octant 1$ and to prove 
the continuity of 7^ -1 on B n D fl Iq. 

For Z = re , the point projected onto the border of E is 


/ __ y _- Y \ 

( rV 1 — e 2 cos-—, rV 1 — e 2 sin-—, er ) . (37) 

\ r cos - r cos - / 

On the other hand, let (X,Y,er) <E P+. From (J26]), (f27|) and ()25jl . its image 
on the sphere, T n (X, Y, er), will be 


x = r\J 1 — e 2 
y = r\J 1 — e 2 


Y 7T 

COS-7-T-, 

n sin - • (X cos - + Y sin -) 
Fvr 

sin- t-^ 5 

n sin - • (X cos - + Y sin -) 

n V n n) 


z = er. 
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But ( X , Y, er) is also a point of ® n , so it satisfies (12U1) . and finally, replacing 
R n from m, the above formulas transform into 


x = r\J 1 — 

y = r\J 1 - 
z = er, 


Y 


£ z COS 


r cos ■ 


Y 


£ z sm 


7 r 5 

r cos - 

n 


i.e. the same as in formula (137|) . 

Similarly one can prove the same equality for a point in IP+ fl® n and then 
in fl B n . In conclusion, the map T~ x is continuous, therefore T n is also 
continuous. ■ 


3.5. Uniform, grids on the sphere 

The most important application of our area preserving projection maps 
is the construction of spherical grids by mapping grids from the polyhedrons 
K n (r, e) to the sphere § 2 (r). The major advantage is the possibility of con¬ 
struction of uniform and refinable spherical grids, since uniform and refinable 
grids are easier obtainable on polyhedrons. 

In the particular case n = 4 and £ = 0, the polyhedron K 4 (r, 0) is an 
octahedron, but not a regular one, as in [I|. However, one can easily see that, 
if we make uniform triangular refinements by dividing each planar triangle 
into four small triangles by lines parallel to the edges, one obtains exactly 
the grids in ||4|. 

Next we will show that the HEALPix grids (see pj) can also be obtained 
by mapping onto the sphere grids on K n (e, r). However, one advantage of our 
map is that it allows us to transport onto the sphere some already constructed 
continuous wavelet bases on polyhedrons, by a similar technique as the one 
described in [3], Compared with the construction in (?;!], one can use now the 
usual scalar product in L 2 (§ 2 ) instead of a weighted scalar product, since the 
map is area preserving. Moreover, the transportation of continuous functions 
on K n (£, r) yields continuous non-distorted functions on the sphere, and this 
does not usually happen when functions defined on rectangles are mapped 
onto the sphere with the area preserving maps used in [1]. When transporting 
a function / : [0, a] x [— b , b] —> M onto the sphere using the maps used in Q, 
the transported function defined on § 2 is usually distorted, and it remains 
continuous on § 2 only when it satisfies some periodicity conditions (see [?{]). 
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Figure 4: The grids on the polyhedron Kg(r,p/(p+ 1)), for p = 2 (left) and p = 4 (right). 
The dashed lines represent the edges of the prism. 

3.5.1. Obtaining the value of £ 

In order to obtain the HEALPix grids on the sphere, we consider r = 1 for 
simplicity, and we construct the grids on the polyhedron K n (l,e) as follows. 

Let p G N. For p even, on each face of the prism B n we first draw 
horizontal lines to divide each face into p/ 2 rectangles having the same area. 
Then, for each of the rectangle we draw the diagonals. Together with the 
triangular faces of the pyramids of K n , they will form a grid consisting in 
n ■ p cells of rhombic shape (see Figure [4j), and we will impose that all this 
cells have the same area A ce n, and this area will be the same after projecting 
the grids onto the sphere. Therefore we have 

Tt 

*4.(E e ) = 47re = npA ce u, -4(§+) = 27 t(1 - e) = -Aeii- (38) 

Elimining A ce n from these equalities we obtain e = ^y. For p odd, we 
first rotate the prism P“ with 7i/n, around the axis Oz and on each face of 
the prism we first draw (p — l)/2 horizontal lines, such that the resulting 
rectangles have the same area, except the “lowest” ones, which have only 
half of the area. Then we draw the rhombic cells as for the case p even, 
except the “lowest” part of the prism (see Figure [6]). The equalities in ()38l) 
also hold, so in this case we must have again e = p/(p + 1). 

Further, for fc 6 N, we perform a subdivision of each rhombic cell of the 
grid on K n (l ,p/(p + 1)) into k 2 cells with the same area, by drawing k — 1 
equidistant parallel lines to each of the edges (see Figure 3). 
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Figure 5: The grids on the polyhedron K 6 (r,p/(p + 1)), for p = 3 (left) and p = 5 (right). 
The dashed lines represent the edges of the prism. 




Figure 6: The refined grids on the polyhedron K 6 (r,p/(p+l)), obtained after a subdivision 
with k = 2, for p = 2 (left) and p = 3 (right). 
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Next we write the equations of the lines of our grids and we show that 
they coincide with the ones deduced in [2] for the HEALPix grids. 

3.5.2. The grid on the pyramid 

We parametrize the triangle with vertices 

(0,0), (R n cos a*, R n sinaj), (R n cosa i+1 , R n sma i+1 ) 


as 


X = u ■ R n cos aj + v ■ R n cos a i+ i, 

Y = u ■ R n sinaj + v ■ R n sin Q!j +1 , 

with u G [0,1], u£ [0,1 — u\. For the parametrization of the face of the pyra¬ 
mid contained in /+, we need to express Z from the equation of the plane con¬ 
taining the points (0, 0, b n +e), (R„ cos cq, R n sin cq, e), (R n cos a i+1 , R n sin a i+ i,e). 
Thus, we obtain 


Z = e + b n {l — u — v). 


The spherical images have the equations 


x 

V 

z 


Vl — z 1 2 
\J\~ z 2 


cos 

sin 


' 2 vr 

V 

n 

u + v 

2 t r 

V 

n 

u + v 



1 — (1 — s)(u + v) 2 . 


J 


) 


If we map the curves u — Ci, v G [0,1 — u\ and v = C 2 , uG [0,1 — v] (with C 12 
constants chosen properly), we obtain the spherical curves in the HEALPix 
grid. Indeed, if we take for example n — 4 and p = 2 (i.e. £ = 2/3), the 
above equations give 


z 

6 


1 - ^(w + u) 2 , 

71 V 71 i 

2 ' u + v + ~2’ 


(39) 

(40) 


where 6 = arctan (y/x). Let us consider 9 t — 6 mod f = f • From (j39j) 
and dUD we obtain 


2 = 1 



7T 

26 t -n 


2 


? 


and also z 
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Now, if we consider u = | =: uz, for each £ G {0,1,..., k} we obtain the 
family of curves 



and this is exactly the equation (20) from [ 2 J. Note that N si ^ e and k in [ 2 ) 
are in our notations k and £, respectively. 

The other class of curves of our grid are 


£ 

V = k =:V ‘ 



which is exacly the equation (19) from [2], with V£ = k/N s ^ e . 


3.5.3. The grid on the prism 

The point C/ePJ’O B n D If of coordinates 


x v 


Vu 


£\ 2m £ 2(i + l)7r 

- R n cos-h —R n cos- 

k J n k n 



. 2m £ . 

R n sm-b — R. n sm 

n k 


2 {i + 1 ) 7 r 

* 

n 


Zu = e, 


and Q 6 P n fl 1„ fl I t of coordinates 

X Q 

VQ 
Z Q 

belong to the same line of the grid of the prism. A point ( X , Y, Z ) of the 
grid of the prism satisfies the equation of the line passing through U and Q 

X — xu _ Y — yu _ Z — zu _ 

X Q ~ x u Vq ~ Vu zq ~ Zu 


£ p 


2ni 


= 1 — — =F - -Rn COS-b — ± - ) R n cos 


n 


£ p 


/1 ^ p\ r > ■ 


k 


£ p 

- ± - ) Rn sin 


2(i + 1)7T 
n 

2 (i + 1)7T 
n 
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When we calculate its image on the sphere using formulas (T5Tj) . (132]) and 
m, we obtain 


z = £ — 2 et, 

x = VTZr ** cos {v'{ i+ i ±P i))’ 

y = VT^ sin + 


In the particular case p — 2, n — 4, £ = 2/3 we obtain 


0 


2 

3 

7T 

2 


4t 

3~’ 


i + — ± t 


If we take again 6 t = 9 mod | = | • (| ± t), we further obtain 

'2 9 ± _t_ 

n k 


2 4 

' _ 3 ^ 3 


which agree with formulas (22), (A2), (A3) in [2j. 


4. Volume preserving map from the interior of a polyhedron K. n (r f , e) 
to the ball of radius r. 


We will use the area preserving map constructed in the previous sec¬ 
tion for the construction of a volume preserving map from the interior of a 
polyhedron K n onto the ball of radius r. 

We fix e, and for p > 0 we define 


K n (p, e) = int(IK n (p, e)) U K n (p, e) and 

&(p) = {(*, y, z) G R 3 , x 2 + y 2 + z 2 < p 2 }. 
Their volumes are 

vol(K n (p,£)) = p 3 -/, woZ(§ 2 (p)) = p 3 /3, 


with 


, c{e) \ TV 2 71 
7 = 2(e + — cot-, 

i n n 


c(e) = \/4(l — £) 2 - 


7P 




cot 2 


TT 


n 


P = 


47T 

T‘ 
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We also need to define the domains 


§+0) = § 2 (p) n {{x,y,z), z > ep} : 

E(p) = § 2 (p) n {(x,y,z), \z\ < £p}. 

4-1. Construction of the volume preserving map V n : § 2 (r) —> K n (ry//3/'y, e) 
Next we fix r > 0, n E N, n > 3, and we will investigate the possibil¬ 
ity of construction of a map V n : § 2 (r) —> 'K n {r-^//3/ r y, e) which is volume 
preserving, i.e. 

vol(D) = vol(V n (D )), for all D C W(r). 


We propose the following construction, which has two steps. 

Step 1. We denote £ = \fj3p) and define the map C : § 2 (r) —» § 2 (£r), 


£{x, y,z) = £- 0, y, z), for all (x, y, z) E S 2 (r). 

Step 2. Let (x, y,z)E § 2 (r£). Then (x, y,z) E 8 2 (p), with p = ^/x 2 + y 2 + i 2 
£p. Now we consider the area preserving map T n ,p '■ S 2 (p) —> K n (p, e) defined 
in Section [3] (in formulas (IT9l) - (l2Tl) . fl22|l - f[2^1) and (T36l) we have to take 
p instead of r). 

Finally, we define the map V n : S 2 (r) —> K n (r£//3/7, e) by 


V n (x,y,z) = (T„,pO £)(x,y,z) = {T n ^ x 2 +y 2 +z 2 p /2 o C)(x, y, z). 

Next we investigate wether the Jacobian of V n ca be ±1. 

Case 1. For (x, y, z) E § + (r) D Iq the formulas for (X, Y, Z ) = V n (x, y, z) 
are 


X 

Y 

Z 


e 


p(p - z ) ( 71 


n sm ■ 


tv y 

— sm — arctan — 


n 


x 


C p(p - Z) cosg.arctan^, 
l — e n x 


£p £ + c(e) 1 


1 ~z/p 
1 — £ 


(41) 

(42) 

(43) 


where p 
we find 


\Jx 1 + y 2 + z 2 . If we evaluate the Jacobian of V n 


restricted to Iq 


J{V n ) 


ax ax ax 

dx dy dz 
dY 8Y dY 
dx dy dz 

az^ az_ az_ 

dx dy dz 


c{s) + £ 

(1 — £)(c(£) + 3£) 


> 0 . 
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If we impose the volume preserving condition J(V„) = 1 we obtain the solu¬ 
tion e = 0 for all n > 3. Another case when J(V n ) = 1 is when c(e) = 2 — 3e, 
which reduces to the equation 


7T n 7T 


2 — 3e = y 4(1 — e) 2 -- cot 2 —. 


r n 


(44) 


This equation gives solutions only in the case when n = 3,4, 5, and the 
solutions are 


for n — 3, 


£ = ^18 — ^3(108 - 5tt 2 )) , i.e. e = 0.105226...; 


for n — 4, 


e 1)2 = — ^8 ± v 7 64 - 5tt 2 ) j , i.e. = 0.20861..., e 2 = 0.59139...; 


for n — 5, 


£ 1j2 = — ( 10 ± yOOO - 5 tt 2 - 2 a/57t 2 ) ) , i.e. = 0.297912..., £ 2 = 0.502088. 


Case 2. For (x, y, z ) G E(r) D 1$ the formulas for (A", Y, Z ) = V n (x, y, z ) are 


A" = £p (-— — sin — arctan — 

V n sin - n x 

\ n 

Y = ip cos — arctan —, 
n x 


z = iz. 

In this case we calculate that 


J{V n ) = 


3e + c(e) ’ 

and its value is 1 when c(e) = 2 — 3e. 

Of course, it is immediate that the Jacobian is 1 for the whole ball § 2 (r). 
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Remark 5. The method that we have used above, which makes use of the 
area preserving map, can be applied only in the cases when all the polyhedrons 
K n (p, e), with p G (0, r£] admit inside a sphere which is tangent to all the 
faces. Indeed, if we split ~K n (p, e) with e 0 into small pyramids with apex 
at the origin O and bases the faces of K n (p, e), there exists such a sphere if 
the distance from O to a face of IP+ equals to the distance from O to a face 
of B„, i.e. 


{b n + ep)Tl 71 Tip 7T 

—-— cot - = — cot - =: d n . 

2(1 — £)n n n n 

In this case, the sum of the volumes of the small pyramids is 


(45) 


_ u _ d n A(Kn(p, e)) _ d n A(S 2 (p)) _ d n 47rp 2 
vol(lK n [p, ejj g g g 

This equals the volume of the sphere E> 2 (p/f) if d n = p/I, 3 , which reduces to 
c{e) = 2 — 3e. 

On the other hand, condition (1451) further gives b n = p {2 — 3e), i.e. again 
c(e) = 2 — 3e. This condition is exactly the condition obtained by imposing 
J (Vn) = 1- 

We remind also the construction in 1 121 . where a similar method worked 
for the cube, which is also a polyhedron which can be circumscribed to a 
sphere. 


4-2. The inverse V~ l : K n {r-^//3/ 7 , £■) —* § 2 (r) 

4-2.1. The case e = 0 

Let (A", Y, Z ) G 1K n (rf//3/ 7 , 0). First we have to find the domain I/ which 
contains (X,Y,Z). Suppose, for simplicity, that (X,Y,Z) G Iff. We have to 
find p G (0,r£] such that (X,Y,Z) G K n (p, 0). The equation of the plane 
containing the face of P+ D if is 


7r 


X + Y tan —V Z 


7r 


7T p 


n 


n sin - \ A — \ cot 2 - 


nsm 


(46) 


Therefore, when (A", Y, Z) is given, then p can be calculated from formula 
(gsp . Further, Vf 1 {X,Y, Z) is defined as 

Vp(X, Y, Z ) = (£-' o T-d)(X, Y, Z), (47) 

which can be immediately be calculated using the formulas fl25]l . fl26|l . f|27]l for 
a.i = 0,£ = 0, and r =~p. The formulas for the general case (X,Y,Z) G if 
cand be deduced similarly using the same ideas as in Section [3] 
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4-2.2. The case £ ^ 0 

For simplicity, suppose again that (X,Y,Z) E Iq U If. In this case we 
have to find not only the radius p of the polyhedron K n (p, e) which contains 
the point ( X , Y,Z), but we must also find whether ( X , Y, Z) belongs to the 
pyramid P+(p, e) or to the prism ® n (p, e). First we observe that the pyramids 
P+(p, e), p > 0 are situated above the conical surfaces generated by lines 
that make a constant angle /3 = arccose with the axis OZ. This means that 
(. X , Y, Z) E P+(p, e) if it satisfies the equations 

X 2 + Y 2 < f^- 1 ) z2 ’ Z >0. 

and (X, Y, Z) E P“(p, £) if it satisfies the first inequality and Z < 0. If 

X 2 + Y 2 > 1^ Z 2 , 


then (X,Y,Z) E B n (p, e). 

Further we have to calculate p. In the case when (X, Y, Z) E P±(p, e), the 
determination of p was already described for the case £ = 0, but in this case 
p is calculated from the formula 


X 


Y tan — 

n 


7T 


c(e)n sin - 


= P 


c(e) 


+ 1 


7T 


n sm - 


So it remains to consider the case (X,Y,Z) E ® n (p,e). The equation of 
the vertical plane which contains ( X, Y, Z) can be written as 


X + Y tan — = —-r—y, (48) 

n n sm - 

n 

and this formula allows us to calculate the value of p. 

Finally, after the determination of p, the inverse can be calculated with 
formula (l47lb 

4-3. Transporting uniform and reftnable grids from the polyhedrons K n (r', e) 
onto the ball 

The most important application of the volume preserving map V n is that 
it allows the construction of uniform and refinable grids on the ball starting 
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from similar grids of the polyhedron. In a first step we construct 2 n tetra¬ 
hedrons (triangular prisms), each having the vertex O and base a face of the 
pyramids P±. For the prism B„, we triangularize each face of B n and then 
construct again 2 n tetrahedrons with vertex O and base a triangle. In the 
second step, we split every of the 4n tetrahedrons into four smaller tetrahe¬ 
drons having the same volume, this simple procedure being described in Q. 
The refinement can be repeated, therefore we end up with a uniform and 
rehnable grid of K n (r', e), whose image on the ball S 2 (r) will be a uniform 
and rehnable grid. 
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